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ABSTRACT 

Aims. We have studied the bulge and the disk kinematics of the giant low surface brightness galaxy ESO 323-G064 in order to 
investigate its dynamical properties and the radial mass profile of the dark matter (DM) halo. 

Methods. We observed the galaxy with integral field spectroscopy (VUT/VIMOS, in IFU configuration), measured the positions of 
the ionized gas by fitting Gaussian functions to the [O III] AA4959, 5007 and 11/3 emission lines, and fit stellar templates to the galaxy 
spectra to determine velocity and velocity dispersions. We modeled the stellar kinematics in the bulge with spherical isotropic Jeans 
models and explored the implications of self consistent and dark matter scenarios for NEW and pseudo isothermal halos. 
Results. In the bulge-dominated region, r < 5", the emission lines show multi-peaked profiles. The disk dominated region of the 
galaxy, 13" < r < 30", exhibits regular rotation, with a flat rotation curve that reaches 248 ± 6 km s"'. Erom this we estimate the total 
barionic mass to be Mi,ar ~ 1.9 • 10" M© and the total DM halo mass to be Mom ~ 4.8 • 10'^ Mq. The stellar velocity and velocity 
dispersion have been measured only in the innermost » 5" of the bulge, and reveal a regular rotation with an observed amplitude of 
140 km s"' and a central dispersion of cr = 180 km s"' . Our simple Jeans modeling shows that dark matter is needed in the central 5" 
to explain the kinematics of the bulge, for which we estimate a mass of » (7 ± 3) ■ lO'" M©. However, we are not able to disentangle 
different DM scenarios. The computed central mass density of the bulge of ESO 323-G064 resembles the central mass density of 
some high surface brightness galaxies, rather than that of low surface brightness galaxies. 

Key words, galaxies: kinematics and dynamics - galaxies: spirals - galaxies: individual: ESO 323-G064 



1. Introduction 

Galaxy properties define a continuum, in size, luminos- 
ity, mass, and surface brightness. Low surface bright- 
ness (LSB) galaxies have significantly extended the range 
in surface brightness over which galaxies can be stud- 
ied, and thus they have received a great deal of atten- 



tion in studies of dark matte r (e.g. [P fenniger et al. 1994 



de Blo k & McGaugh 1997; McGaugh et al. 2001; Swate rs et al. 
2000, 2003; Kuzio de Narav et al. 2008), stellar content (e.^ . 
Ide Blok et al., ,1995i: ,Bothun et al.. ,1997; 'Boissier et al.' '200? 



and ga s content (e.g. lO'Neil & Schinnere r 2003; Pizzella et al] 
l2008 b). Because of their unique properties, LSB galaxies also 
play a significant role in our understanding of the universe. Their 
distribution in the field and in galaxy clusters will help us to 
understand the bright and dark matter distribution in the lo- 
cal neighborhood and in the universe; their number and density 
distribution may be relevant to the stud y of damped Lv-g ab- 
sorpti on seen against quasars (e.g. .Jimenez et al.lll999t O'Neill 
I2OO2I) . Like high surface brightness galaxies, LSB galaxies also 
span a range in properties, such as size, mass, and bulge size. 



Send offprint requests to: L. Coccato, 
e-mail: lcoccato@mpe.mpg.de 

* Based on observations carried out at the European Southern 
Observatory (ESO 075.B-0695). 



For example, 'Bothun et al.' (I99(f} discovered a class of giant 
LSB galaxies, and Beiiersbergen et al. ( 1999a) studied a sample 
of bulge-dominated LSB galaxies. To date, very little is known 
of the properties of giant-sized or bulge-dominated LSB galax- 
ies even though they can enhance our understanding of galaxy 
formation and evolution by providing an opportunity to sample 
galaxies in a hitherto almost unexplored regime of galaxy prop- 
erties. Several scenarios have b een proposed to expla in their for- 
mation: density peak in voids ("Hoffman et al."1992'), bar insta- 
bility (|Npguchi 2001; Maver & Wadslev 2004), or secular evo- 
lution from ring galaxies (IMapelli et al.ll2008n . 

To fill the gap between regular and giant LSBs, Swaters, 
Rubin, & McGaugh (SRM in preparation) have recently com- 
pleted a study of the kinematics of a sample of bulge dom- 
inated LSB galaxies. Their aim is to study the dark mat- 
ter properties and to determine whether bulge-dominated LSB 
galaxies are dark matter dominated like other LSB galaxies 
(Ide Blok & McGaugh|[T997tlSwaters et ai]l200l . 

In this paper we show additional results obtained for the gi- 
ant LSB galaxy ESO 323-G064, selected from the SRM project. 
ESO 323-G064 has an heliocentric velocity of 14830 km s^' at 
a distance of 194 Mpc (assuming Hq = 75 km s"' Mpc"' and 
a galactocentric correction of -278 km s"'). On its digital sky 
survey and 2MASS images, ESO 323-G064 has a bright, com- 
pact bulge, suiTounded by a low surface brightness disk. There 
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is some evidence of a weak bar component. The long slit spectra 
of ESO 323-G064 obtained (SRM, in preparation) with the 6.5m 
Baade Telescope at LCO, show strong, double-peaked emission 
from Ha, [N II], [SII], [OI] in the nuclear region. The large ex- 
tent (47 kpc) of the Ha and the large peak-to-peak rotation veloc- 
ity of 445 km s^' (on the plane of the sky) mark ESO 323-G064 
as a likely giant LSB. 

In order to determine the kinematics of the stars in the bulge- 
dominated region, and to compare these motions with those of 
the gas in both the nucleus and at larger radii, we decided to 
obtain two dimensional data using an integral field unit. This 
will give a better picture of the complex gaseous kinematics 
especially the double peaked emission revealed by the long 
slit data. Optical two-dimensional observ ations in LSB galax- 
ies have been ob tained in the past (i.e. 'Swaters et al.J 20031: 
iKuzio de Naray et al. 2006, 2008; Pizzella et al. 2008b) but, so 
far, never for the stellar component of a giant LSB. 

The paper is organized as follow: in Section|2]we present the 
integral field observations and discuss the data reduction pro- 
cess; in Sections 3 and 4 we present the kinematical results for 
the ionized gas and stellar components respectevely; in Section 
5 we present the Jeans model to the stellar kinematics and in 
Section 6 we discuss the results. 

2. Observations and data reduction 

The integral-field spectroscopic observations were carried out 
in service mode with the Very Large Telescope (VLT) at the 
European Southern Observatory (ESO) in Paranal (Chile) from 
June to August 2005 during dark time. The Unit Telescope 
3 (Melipal) was equipped with the Visible Multi Object 
Spectrograph (VIMOS) in the Integral Field Unit (IFU) config- 
uration. The seeing, measured by the ESO Differential Image 
Meteo Monitor, was generally below (J.'! except for observa- 
tions on July during which the seeing was around r.'5 (one 
30 minute exposure) and observations on 4* August in which 
the seeing was between 0'.'7 and I'.'O (one 30 minute exposure). 

The observations were organized into 8 exposures of 30 min- 
utes each, divided into two different pointings with an offset of 
25 arcsec. Each telescope pointing has a field of view of 27 x 27 
arcsec^ and it is recordered on a 4 CCD mosaic. We hereafter 
refer to the 4 CCDs as quadrants, and we refer to the 2 observed 
fields as field A and field B. Quadrants # 1,2,3 and 4 cover the 
NE,SE,SW and NW portions of the field of view respectively. 
Field A (4 X 30 min exposures, with a few pixel dithering) covers 
the NW side of the galaxy, while Field B (4 x 30 min exposures 
with a few pixel dithering) covers the SE side of the galaxy. The 
field of view of the four VIMOS quadrants was projected onto a 
micro lens array. This was coupled to optical fibers which were 
rearranged on a linear set of micro lenses to produce an entrance 
pseudo-slit to the spectrograph. The pseudo-slit was 0'.'95 wide 
and generated a total of 1600 spectra covering the field of view 
with a spatial resolution of 0'.'67 per fiber Each quadrant was 
equipped with the HR-blue resolution grism (4120 - 6210 A) 
and a thinned, back-illuminated EEV44 CCD with 2048 x 4096 
pixels of 15 X 15 jum^. The spectral resolution measured on the 
sky emission lines is A/AA » 1700; cr,„j„. x 75 km s"'. 

Together with every exposure a set of night calibration spec- 
tra were taken: one comparison spectrum (Neon plus Argon) for 
the wavelength calibration and 3 quartz lamp exposures for the 
flat field correction and fibers identification. 

For each VIMOS quadrant all the spectra were traced, iden- 
tified, bias subtracted, flat field corrected, corrected for rela- 
tive fiber transmission, and wavelength calibrated using the rou- 



tines of the ESO Recipe Execution pipeline (ES0rex)[3. Cosmic 
rays and bad pixels were identified and cleaned using standard 
MIDAS ' routines. We checked that the wavelength rebinning 
was done properly by measuring the difference between the 
measured and predicted wavelength for the b rightest night-sky 
emission lines in the observed spectral ranges jOsterbrock et al.l 
11996) . The resulting accuracy in the wavelength calibration is 
better than 5 km s The intensity of the night-sky emission 
lines was used to correct for the different relative transmission 
of the VIMOS quadrants. The processed spectra were organized 
in a data cube using the tabulated correspondence between each 
fiber and its position in the field of view. From the single expo- 
sures of every observed field we built a single data cube. The 
spectra were co-added after correcting for the position offset. 
The offset was determined by comparing the position of the in- 
tensity peaks of the two reconstructed images obtained by col- 
lapsing the data cubes along the wavelength direction. The accu- 
racy of the offset is +0.5 pixel (^ 0'.'33). This slight deterioration 
of the spatial resolution does not affect the results. 

Finally, we co-added the two available data cubes (field A 
and field B) by using the intensity peaks of the flux maps as 
a reference for the alignment. In this way we produced a single 
data cube to be analyzed in order to derive the surface brightness 
and two dimensional field kinematics. 

In Figure [T] we show the two reconstructed images of both 
the observed fields, obtained by collapsing the data cube along 
the direction of dispersion. 

2.1. Sky subtraction 

The sky subtraction is a critical step in the data analysis (es- 
pecially for stellar kinematics) because 1) the galaxy absorption 
lines are weak and 2) some weak sky emission lines overlap with 
the galaxy emission lines. 

The automatic ESOrex pipeline evaluated and corrected the 
spectra for the different efficiency of the fibers and the 4 de- 
tectors. But in order to avoid possible problems due to a poor 
correction of the efficiency difference for the four quadrants, we 
decided to evaluate the sky contribution for every quadrant sepa- 
rately. This choice is supported by the fact that the useful signal 
for the stellar kinematics is contained entirely in one quadrant 
(quadrant #3 for field A observations and quadrant #4 for field B 
observations, respectively). 

For each quadrant we identified the spectra in which the 
galaxy's contribution is negligible and we computed the median. 
Finally we subtracted the median sky spectrum from every sin- 
gle spectrum of the quadrant. 

3. Gaseous kinematics 

ESO 323-G064 presents very bright emission lines in the nu- 
cleus together with fainter emission lines located at larger radii. 

The gas kinematics are measured by fitting to each spectrum 
a background level and several Gaussian components. The num- 
ber of components is chosen depending on which and how many 
emission lines are visible at the considered spectrum. No con- 
straints are applied to the fitting parameters: velocities, velocity 
dispersions and intensities are independent. 

To determine the kinematics, we used a non-linear least- 
squares minimization algorithm bas ed on the robust Le venberg- 
Marquardt method implemented bv lMore et alj ( Il980l) . The ac- 

' ESOrex and MIDAS are developed and maintained by the European 
Southern Observatory 
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Fig. 1. Left panel: R-band image of ESO 323-G064 from the ESO-Uppsala Galaxy surface photometry catalog (ESO-LV. iLaubertsI 
[l982). North is top, east is left. Right panel: reconstructed image of ESO 323-G064 from combination of 8 exposures. Scale and 
orientation are as in the left panel. 



300 - ^ 
250 |f^^* 




Fig. 2. A: Rotation curve of ESO 323-G064 from HyS emission, formed from the two dimensional velocity field of the disk (open 
diamonds). Filled circles are the stellar rotation circular velocity curve after correction for asymmetric drift (see 15.4. It . B: Velocity 
field of HyS emission at large nuclear distance. C: Velocity dispersion field of HyS emission at large nuclear distance. In both panels 
B and C, galaxy central contours are shown, the field of view is 46" x 46", North is at top. East at left. 



tual computation has been done using the MPFIT procedure im- 
plemented by C. B. Markwardt under the IDL environment 



3.1. Ionized gas at large radii 

Faint HyS emission is detected at 15" < r < 35". At this distance 
from the center, the gas is associated with the disk component. 
It reaches ^ 15050 km s"' on the NE side and ^ 14500 km s"' 
on the SW side. The intrinsic velocity dispersion ranges between 
30 and 50 km s ' showing that this gaseous component is quite 
cold (see Figure|2|. 

To derive a rotation curve from the two-dimensional velocity 
field of the gaseous disk, we divided the gaseous disk in concen- 



^ The updated version of this code 
|http://cow.physics.wisc.edu/~craigm/idl/idl.html| 



is available on 



trie rings, with the same inclination (/) and the same orientation 
on the sky (PA). Each ring is constructed in order to contain the 
same numbers (20) of data points, with velocities v„ and posi- 
tion angles (f>„ B We also assume the gas is moving in circular 
motion, therefore the velocity VMiR) in the M-th bin centered 
at R is measured by fitting the following function to the N - 20 
points within the bin: 



V((p,R) = Vm{R) ■ cos(0 + PA) ■ sin(0 + V„ 



(1) 



The best fit results in a systemic velocity of Vjvj = 14830+12 
km s"', a position angle on the sky of PA = 38 + 3 (slightly dif- 
ferent from the 32 degrees derived on the ESO-LV R-band and 
B-band images, an inclination of / = 61.5 ± 0.9 (consistent with 



p„ are measured counterclockwise starting from the galaxy major 
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the 60 degrees from ESO-LV images Q.) and a center consistent 
within a pixel with the photometric center of the bulge. At ev- 
ery bin position the IDL fitting procedure gives the amplitude 
of rotation Vm{R) and its error. We show in the left panel of 
Figure|2]the derived rotation curve of ESO 323-G064. The rota- 
tion curve is relatively constant with radius. The average rotation 
value (taking into account the incliation) is 248 ± 6 km s"'. The 
result does not change if we hold the position angle constant at 
PA - 32 as derived from ESO-LV images. 

Since the velocity dispersion of the ionized gas is quite small, 
we can approximate the measured velocity rotation with the cir- 
cular velocity Vc- This allows us to estimate the total baryonic 
mass of the galaxy, from the baryonic TuUy Fisher relation, as 
described in lMcGaughl (|2005|) : 

Mbar = 50. ■ V^. = (1.9 ± 0.2) ■ 10"Mo (2) 

In addition, we can also make a rough estimate of the total 
mass Mom of the dark matter halo through the relation: 

Mdm = 2.52 ■ 10'- (^)' = (4.8 ± 0.4) ■ lO'^M^ (3) 

To derive the formu la, we followed the pre scriptions by 
iBrvan & NormanI (Il998l) and lBuUock et alj (1200 1[) using Qa = 
0.7, Q„ = 0.3, h = 0.7, A, = 105.4 and z 0.05. The error 
estimates for Mt,ar and Mdm are simply given by the error prop- 
agation in Equations |2] and [3] 

3.2. Ionized gas in nuclear region 

Hy8, [O III] /14959 and [O III] /15007 bright emission lines are vis- 
ible in the nuclear regions within r < 5". Interestingly, at some 
locations these emission lines are triple peaked (see Figures [3] 
and|4|. This complexity is probably generated by three barely 
resolved emission-line regions located close to the center but at 
different azimuthal angles. 

In order to have kinematical informations for all the 3 re- 
gions, we fitted each of the the galaxy emission lines with three 
independent Gaussian functions, and the results (velocity, veloc- 
ity dispersion and intensity) were compiled separately. 

The velocity field of the three components do not show a reg- 
ular rotation. Their averaged radial velocities are approximately 
14539 + 38, 14726 ± 31, and 14992 ± 26 km s"'. These values 
were inferred from the [O III] velocity field (mean value between 
[O III] .14959 and [O III] /15007) because it is not affected by stel- 
lar absorption as is the li/3 emission line. The velocity dispersion 
in each of the 3 regions is between 50 and 100 km s"' . 

The measured FWHM of the 3 regions is about 1'.'7, which is 
very similar to the one of the foreground star 1'.'3). The max- 
imum distance measured between the centroids of the 3 regions 
is I'.'O. 

We can exclude that the observed multi-peak line profile is 
an artifact caused by the data reduction because: 

1. The sky lines observed in each quadrant (either in the raw 
spectra or in the final reduced data cube) are not multi 
peaked. This eliminates instrumental effects (e.g., a shift of 
the instrument during the exposure) or pipeline errors (e.g., 
a shift in the dispersion direction during the data cube com- 
bining) as a source of the triple-peaked profiles. 

Images from scarmed plates can be downloaded from 
|http://www. astro- wise, org/portal/aw-datasources . shtml [ 




5180 5200 5220 5240 5260 5280 5500 
wavelength [angstrom] 

Fig. 3. Examples of [O III] emission line profiles, in which the 
3 peaks are visible. H/3 emission line has been omitted in this 
plot because it is very faint with respect to the [O III]. The blue, 
green and red lines represent the fits to the first, second and third 
emitting regions respectively. 

2. The three peaks in the galaxy emission line are present be- 
fore and after the sky subtraction and they are present also 
in the raw spectra. Moreover, the few sky lines which are 
overlap with the galaxy emission lines are very weak. 

3. The three peaks are visible in all the galaxy emission lines 
and the ratio between the intensities of the components 
is (almost) the same if measured in H/3, [O III] /i4959 or 
[O III] ^5007. 

4. The peaks are visible also if we consider only the exposure 
with the best seeing condition. 



4. Stellar kinematics 

The stellar kinematics have been measured by means 
of the Penali zed PixeLFitting (ppxf) method by 
ICappellari & Emse llem (20 04|). We c hose a library of stel- 
lar templates from iValdes et al.l (|2004|) provided together with 
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Fig. 4. Upper panel: position velocity diagram towards the 
galactic nucleus extracted along the major axis. Lower panel: 
intensity profile of the PV diagram along the dashed line in the 
upper panel. The 3 peaks of the emission line are clearly visible. 



the ppfx tools. These spectra have a spectral range of 4780 - 
5460 A and spectral resolution of 1.0 A at 5100 A (tr ~ 22 
km s"', which had been deteriorated in order to match the 
instrumental one). 

The stellar templates from this library were choosen among 
type K stars, which are commonly used as kinematical tem- 
plates. 

Errors of the measurements ranged around 15-20 km s"'. 
We did several tests using additional t emplates from the same 
library and from the MILES library jSanchez-Blazquez et al.l 
|2006|) . Results differ from the original values with a scatter of 
10 km s which is consistent with the measurement errors. 

No multiple peaked features in the absorption line spectra 
are visible at our spectral resolution (and at our signal to noise 
ratio). Some examples of the fit results are given in Figure|5]and 
the two-dimensional velocity and velocity dispersion fields are 
shown in Figure|6] The velocity of the stars ranges from ~ 14700 
to ~ 14900 km s"', and a clear pattern of rotation is visible. 

The central intrinsic velocity dispersion within 1" is a; 180 + 
10 km s"', which, combined with the gaseous circular velocity 
value Vc - 248 km s"' (derived in Section lTTT i. set this galaxy 
close the {Vq - ctc) relation for LSB galaxies (Courteau et aQ 
l2007 a b). According to these authors, the predicted Vc value for 
a bulge with cr^ = 180 km/sec is Vq ~ 290 km s"' . 



The Vy a of the bulg e is computed u sing the formahsm intro - 
duced bv lBinnevI (l2005l) and apphed bv lCappellari et all ( l2007h . 



\criiNTR \cr/oBS 



^n=\ « ^ n 



0.56 + 0.02 



(4) 



where F„, y„ and cr„ are the flux (extrapolated from the sur- 
face brightness profile), velocity and velocity dispersion mea- 
sured from the n - th spectra. An inclination of 62 degrees and 
isotropy assumption are adopted to evaluate the measurements 
in the edge-on case. Errors are calculated using classic formulas 
for error propagation. 

This value is consistent with the predicted value of an edge- 
on isotropic oblate system (Binnev 1978) 0.4 < (^)^^^^ < 0.6 
for an intrinsic ellipticity ranging 0.13 < ejNTR < 0.27. The in- 
trinsic value for ellipticity has been calculated from the observed 
one ( 0.1 < eoBS < 0-2, Se ction|5.11l after correction f or inclina- 
tion using the formula (see lBinnev & Tremaine|[T987l) : 



eiNTR 



1 - V^'T6o^^(eo^7^^2j7sin(0 



(5) 



5. Stellar dynamical modeling 



In this section we present a simple dynamical model of the stellar 
velocity and velocity dispersion field^ . 

We assume the total luminous plus dark matter mass distri- 
bution of the central galaxy bulge is spherical. The low reported 
value for the ellipticity (Section [5. II ) allows us to make this as- 
sumption. In addition we assume isotropy for the bulge, given 
the fact that the inferred V/cr is consistent with an isotropic ro- 
tator, as shown in SectionH] 

The main deviation from spherical symmetry in the bulge of 
ESO 323-G064 is the bar component, clearly visible in Figure[T] 
Given the limited amount of stellar kinematic data in the bulge, 
it is not possible to take into account the complexity of stellar 
motions in the bar potential. However, from the isophotes shown 
in Figure|6]it is evident that the bar contribution dominates only 
the outer regions of the velocity field, for R > 3'.'5. Stellar or- 
bits in those regions will be elongated towards the bar direction, 
which is nearly orthogonal to the kinematic major axis. This will 
lead to an underestimation of the stellar rotation velocity and it 
will be quantified in Section|53] 

We assume also that the disk contribution to the total mass in 
the bulge region (inner 5") is negligible. We do not have direct 
measurement of the stellar disk kinematics, but its low luminos- 
ity (compared to the bulge) supports our assumption. Moreover, 
several studies demonstrated that the disk contribution to the in- 
ner dynam ics in low sur face brightness galaxies is almost in- 
significant dSwaters et alJlIOO O. 2003). 

The limited spatial extent of the observations and the lack of 
good photometry for this galaxy do not allow us to go in more 
detail. 

The Jeans equation in radial coordinat es for such a sys tem 
is (we followed the formahsm adopted bv iHui et al.lll995] and 
iPeng et al.ll2004h : 



Pk{r)crfjr) 



GM(r)pk(r) ^ pk{r)Vl,{r) 



(6) 



We assume Ho = 75 km s ' Mpc ' in the modeling. 
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Fig. 5. Examples of the stellar kinematics fit quality. Black is the galaxy spectrum, while red (green) represents the portion of the 
stellar template which was included (excluded) from the fit. 



which, when solved for the intrinsic velocity dispersion cr;„,, is: 

'(x)-xVUx) 
) 7 dx (7) 



1 

= —- Pkix. 

Pkin Jr 



where is the mass density radial profile for the tracer of the 
potential (i.e. the stars), Vroi is the intrinsic rotation curve of the 
stellar component, M is the total mass of the galaxy and G is the 
gravitational constant. 

The velocity dispersion projected on the sky is: 



cT^(R)^Vlo,(R)-VARf (8) 

where is the projected rotation curve of the stars (see Section 
I5.2l i and Vws is the line-of-sight second velocity moment, pro- 
jected into the sky, given by: 



2R 



Pk(r)r 



r-l Vr2 - 



.dr 



(9) 



5.1. Density of the kinematical tracer 

There are no good CCD images of ESO 323-G064 available in 
the literature. The better available photometry is from the ESO- 
LV catalog. 3 We perform ed a bulge/disk decomposition using 
galfit tPeng et al.ll2002h . adopting a de Vaucouleurs law for 



the bulge and an exponential law for the disk. We derived an 
effective radius R^ - Ql.'bl and an average ellipticity of 0.3 for the 
bulge, position angle PA - 32 and inclination / - 60 for the disk. 
With the ESO-LV images we confirmed also that the stellar disk 
of ESO 323-G064 is in the LSB regime. Its surface brightness 
in the B-band ranges from ~23.3 mag arcsec"^ at the center to 
~26.5 mag arcsec"^ at 35". Even though the ESO-LV images 
have been useful to get an estimate of the surface brightness, 
they had insufficient signal-to-noise and spatial resolution to get 
reliable values of the bulge parameters. We therefore used our 
VIMOS observations, collapsing the observed data cube along 
the dispersion direction and deriving a surface brightness profile 
from the resulting image. We fit a de Vaucouleurs law to this 
profile and derived an effective radius Rg - 0'.'9 (see Figure |7]) 
and a ellipticity ranging from 0.1 to 0.2, not too far from the 
adopted spherical approximation. Although the signal-to-noise 
ratio of the VIMOS observations is low, the advantage is that the 
VIMOS spatial resolution (0'.'67 / pixel) is higher than that of the 
ESO-LV images (r.'35 / pixel). This makes us more confident 
of bulge parameter values based on VIMOS observations than 

values based on ESO-LV images. 

The de Vaucouleurs profile resembles the iHernquisl (Il990h 
mass model in which the intrinsic density distribution is given 
by: 



Pk(r) = 



1 



27r rir + df' 



(10) 



^ See Section lSTTI for the images source. 



where Mi is the total luminous mass and a - /?(./L8153. 
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Fig. 6. Stellar radial velocity and velocity dispersion fields, with 
the bulge central isophotes shown. The field of view is 14" x 14", 
North is up. East is left. 



The projection on the sky of Equation [TO] is given by: 



f 

Jr 



Pk(r)r 
Vr2 - 



(11) 



Equations[TO]and[TT]will be used in Equation|9]to determine 
the velocity moments for the mass calculation. 




LO l^^ 1,4 

Ri/4 (arcseci/4) 



L6 



Fig. 7. Dots: Surface brightness profile obtained by collapsing 
the data cube along dispersion direction. Continuous line: de 
Vaucouleurs fit to the measured surface brightness, with = 
0'.'9 and /le - -16.8 mag arcsec"^ (arbitrary zero point). The in- 
nermost point was not included in the fit, since it is affected by 
the seeing. 

5.2. Intrinsic and projected rotation curve 

Because of the limited spatial extent and resolution of the obser- 
vations, a "bidimensional" approach to the Jeans equation can- 
not be reliably carried out. We therefore decided to average az- 
imuthally and then bin radially the bidimensional velocity and 
velocity dispersion fields. To do that, we divided them in ra- 
dial bin^ (after several tests, it turned out that I'.'S per bin is an 
optimal value) and for each radial bin at a given position we: 

- calculated the velocity Vir^) by fitting the function v'(0) = 
Vir^) ■ cos(0 + P.A.) + Vsys to all the velocity data points 
within the N — th bin. is the angle measured in the galaxy 
equatorial plane, P.A. is the galaxy position angle on the sky 
and Vsys - 14825 km s"' is the adopted systemic velocity. 

- calculated the velocity dispersion a-f^ir^) by averaging all 
the velocity dispersion data points within the - th bin. In 
this computation, we took into account also the value of the 
velocity gradient along the resolution element (I'.'S) and re- 
moved it from the velocity dispersion. 

In Equation]?] the quantity V,,,, is the intrinsic stellar rotation 
curve. It can be parametrized with the expression: 



VroAr) = 



4 



(12) 



where is its asymptotic value for the velocity and r/, is a scale 
parameter The adopted function was chosen in order to properly 
match the observed velocity curve with the minimum number of 
free parameters. This expression has been used for example also 
in the spherical PNe system of Centaurus A jHui et al.l 119951 : 
iPeng et al.ll2004. 

The observe d rotation curve is obtained projecting Equation 
[121 on to the sky dBinnev & Tremainelll987l) : 



VsiR) = 



2R 
W) 



f 

Jr 



Pk(r)Vrot(r) 



■R2 



dr 



(13) 



0. 15 in order to take into 



' Bins have an elliptical shape with < 6 >= 
account the bulge ellipticity we observed. 

^ We fixed the position angle to the value P.A. = 38 derived in Section 
13.11 We tested this assumption leaving P.A. free to vary in the fit and we 
found the value to be consistent with 38°. 
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Equations [12] and [13] will be used in Equations [8] and [9] to 



compute the velocity moments. 




150 




5.3. Radial mass profile 


o 

0) 

m 


100 




The radial distribution of the total mass M{r) of the galaxy bulge 
is given by the sum of the luminous and dark matter contents [_]. 
The contribution of the luminous component is obtained by 


V [km 


50 



/ 



integrating Equation[TO]in the volume. It leads to: 



Mkir) 



(r + af 



(14) 



Together with the self consistent case, in which the total mass 
of the galaxy is given only by the contribution of the stars, we 
explored also t wo di fferent scenarios for the dark matter content: 
iNavarro et al.l d 19971 NEW hereafter), and the pseudo isothermal 
halos. 

Their mass density distribution are: 



PNFw(r) = 
Pisot(r) 



Po 



1 + {rlraf 



(15) 
(16) 



where rj, ps, ro and po are scale parameters. 

The corresponding mass distributions are given by the vol- 
ume integration of Equations [15] and [16] We obtain: 



ln(l+r/r,)- 



1 + rlvs 



MisoAr) = ^Tiporl [r/ro - arctan (r/ro)] 



(17) 



(18) 



In the dark matter scenarios, the total mass M(r) is given by 
adding to Mk one of the contributions expressed in Equations]!?] 
or[T8] 



5.4. Model results 

We fit Equations [8] and [13] separately to the velocity Vir^) and 
velocity dispersion o-(rN) data points (see Section [STSI i. The ac- 
tual fit computation was done using the MPFIT algorithm, as 
done in Section[3] 

One could add constraints to the halo parameter by including 
the disk H/3 rotation curve and the disk surface brightness profile 
into the fit. However, given the large uncertainties on the disk 
photometric profiles, we decided to focus our mass model only 
to the bulge regions. For completeness, we present the models in 
which the disk data are used in Appendix A. 

The fit results are shown in Figure [8] and are listed in Table 
[T] Figure [8] shows the best fit model results (solid line) together 
with the observed radial profiles of velocity and velocity dis- 
persion derived in Section 15.21 The last observed point of the 
rotation curve is below the best fit rotation curve, possibly the 
result of the bar influence on the stellar kinematics for R > 
3'.'5 as discussed in Section JS] Assuming this deviation is en- 
tirely due to the bar, a rough estimate of its effect on our results 
can be obtained from a fit of the velocity curve without taking 
into account the last velocity measurement. This is shown as the 



' In Section[5]we assumed that the disk contribution to the bulge mass 
and dynamics is neglegible. 



R [arcHec] 



200 
150 
100 

50 





1 2 3 4 5 

R [arcsec] 

Fig. 8. Upper panel: stellar measured velocity curve of the bulge 
(filled circles, see Section I5.2l i and the best fit model from 
Equation [13] The black line represents the best rotation curve 
model, while the gray line represents the best rotation curve 
model excluding the last data point from the fit. Lower panel: 
measured radial velocity dispersion (filled circles, see Section 
15.21 ) compared to the NFW (solid line), the pseudo isothermal 
(dashed line) and the self consistent (dotted line) best fit models. 



gray line in the upper panel of Figure [8] The maximum differ- 
ence between the two rotation curves (within 4'.'5) is » 14%. 
Mass scales as the square of the velocity, therefore we would 
expect a maximal underestimation of the total bulge mass of 
AM/M = 2Ay/y = 28%, which less than a factor of 1.3. 

From Figure[8]we can see that the self consistent model pro- 
vides a poor fit to the data, compared to the dark matter scenar- 
ios. The bulge kinematics are therefore better explained, under 
our assumptions, with the presence of DM, although we are not 
able to disentangle between NFW and pseudo isothermal mod- 
els. 

The reduced x' of both dark matter models is close to 1 , 
while the self consistent case gives a value of 3.6. In Figures 
12] and [To] we show the Icr, 2cr and 3>cr confidence levels in 
the parameter space. Those plots were produced by scaling the 
measured x^ value to the ideal one x' = N- M- 3(N = 6is 
the number of data points, M = 3 is the number of free param- 
eters) and taking the ex pected A;t^^ variat ions for 3 parameters 
(i.e. 3.53, 8.02 and 14.2: lPress et aL][T99l chap. 15.6). 

The total mass of the bulge (calculated using Equations [17] 
and [18] for r = 5") is (7.4 ± 3.2) ■ 10"' Mq and (7.1 + 3.6) ■ 
10'" Mq according to the NWF and pseudo isothermal scenarios, 
respectively, while the ratio of the dark matter content to the total 
mass is about 0.55 and 0.42, at that radius. In the no dark matter 
scenario the bulge mass is (6.5 ± 1 .6) ■ 10'° Mq. The total fraction 
of dark matter in the bulge, as a function of the radius, is shown 
in Figure [TT] 

Errors of the bulge mass are computed by applying error 
propagation formulas to Equations [T4l[T7] and[r8] 
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Table 1. Bulge parameters from the best fit model. Values are 
converted assuming D - 194 Mpc on the right column. 



Parameter 


Value 




Rotation curve: 


174 ±26 [km s"'] 




rh 


1.6 ± 0.8 [arcsec] 


(1.5 ± 0.8 [kpc]) 


Self consistent 

Mi. 


? = 3.6 
6.6^1;^ [lO'OMo] 




NFW: 

Ml 


? = 0.8 
2.3'4,+ [10'»Mo] 






iOA^l [logio (Mo/"')] 


(p, = 15:^ [Mo pc-']) 


rs 


OJ*qI [arcsec] 


(0.67 ±0.19 [kpc]) 


Isothermal: 

Ml 


? = 1.5 
4.5+3 [IO'^Mq] 




logioPo 


9-6!SJ;i [logio (Mo/"')] 


(po = 5!? [Mo pc-']) 




0.4+j; [ [arcsec] 


(0.38 ± 0.09 [kpc] ) 



5.4.1. Mass density radial profile 

From the measured stellar rotation curve and the radial velocity 
dispersion profiles (i.e. data points in Figure [Hll we derived the 
mass density profiles and compared them with the best model 
predictions. 

First we calculated the circular velocity Vc using the asym- 
metric drift correction (assuming isotropy): 



2 dMy) 



(19) 



If we insert the observed bulge r'^^ radial profile for the light 
distribution v of the kinematic tracer we obtain: 



y2(r^) = vV^) + 1.92(rV)^ 



1/4 



(20) 



Stellar mean circular velocity obtained with the asymmetric 
drift correction is Vc = 243 ± 6 km s"' which is consistent 
with the one measured from the gas in the disk regions (248 ± 6 
km s"'). 

Finally we ca lculated the mass density as done in 
Ide Bloketal.l(l2001h assuming a spherical mass distribution: 



p(rN) 



1 

4^ 



-f +2— 

4 



dVc 
dR 



(21) 



We did not fit these calculated mass density values because 
they do not contain more information than the observed V and 
cr themselves. Moreover, our direct fit to the observed quantities 
does not depend on the assumption adopted in the mass density 
calculation. We show in Figure [T2]the comparison between the 
derived p and the prediction by the NFW and pseudo isothermal 
dark matter halo models. 



6. Discussion and conclusions 

We presented the two-dimensional velocity and velocity disper- 
sion fields for the gaseous and stellar components of the LSB 
galaxy ESO 323-G064. The gas emission lines show a very 
bright and complex structure within the central 5", character- 
ized by 3 peaks, which we interpret as due to the presence of 3 
spatially unresolved emitting regions. At radii out to 30", in the 
region dominated by the galaxy disk, a weak H/? emission is de- 
tected showing a regular velocity field with a maximal amplitude 
of248 + 6kms-i. 

The stellar absorption lines are detectable only in the very 
innermost regions of the bulge. The stellar kinematics shows a 
regular rotation with an amplitude of + 140 km s"' and a central 
velocity dispersion of x ISO km s"'. The intrinsic (V/crjjNTR 
ratio for the bulge is 0.56 ± 0.02, which is consistent with an 
isotropic rotator in the 0.1 < £obs < 0.2 observed ellipticity 
range, considering an inclination of 62 degres. These values of 
y/cr an d € place the bul ge of ESO 323-G064 among fast rotator 
bulges (ICappellari et al.li200 6). The value of V/cr is considerably 
small when compared to the average value determined from the 
sample of 6 bulge dominated LSB galaxies of Pizzella et al. 2008 
(< y/cr » 1). 

The circular velocity, 248 + 6 km s"' measured from the 
gaseous disk, places ESO 323-G064 in good agr eement with the 
locatio n of LSB galaxies in the Vc - cTc plane dCourteau et al.l 
l2007 a''b'). On the other hand, this value is low er then the value of 
320 km s 'predicted by Pizzella et alj (|2005|) . 

The intrinsic bulge ellipticity value for ESO 323-G064 
(eiNTR = 0.20 ± 0.07) is consistent with the mean value of bulges 
of the high sur face brightness dis k galaxies < chsb >- 0.15 (as 
determined bv lMendez-Abreu"era l. 2008) while it is lower than 
the average value for bulge dominated L SB galaxies < eisB >= 
0.45 (as found by iPizzella et"ani2008al in a small sample of 6 
bulge dominated LSB). 

The amplitude of the gaseous rotation curve (248 ±6 km s"') 
leads to an estimate of the total baryonic mass in the galaxy of 
Mbar = (1.9 ± 0.2) ■ 10" Mpusing th e empirical baryonic Tully 
Fisher, as done in llVlcGaughl (120051) . IVIoreover, under the hy- 
pothesis of a ACDM universe (see Section lTTI for details on the 
assumed parameters), we estimate the total mass for the dark 
matter halo Mdm ~ 5 ■ 10'^ Mq. 

We produce spherical isotropic Jeans models for the stel- 
lar kinematics in the bulge, exploring the self consistent, NFW 
and pseudo isothermal scenarios. Even though the data are to 
be taken with some caveats (due to the lack of good photometry, 
limited spatial extension and resolution of the stellar kinematics) 
with this simple analysis we show that dark matter scenarios fit 
the data better than the self consistent model. The derived total 
bulge mass is (7 ± 3) ■ 10'*' but we are not able to disentangle 
between the two different dark matter models. 

The derived central bulge mass density (see Table [U is 
p = 15^5 [Mq pc"^] in the NFW scenario, and p = 5^1 [Mq 
pc~^] in the pseudo isothermal scenario. Typical values of cen- 
to pc" 



tral mass density range from few 10 ^ Mo pc ^ to few 10 ^ 



Ma 



2£. 



^ for regular low surface brightness (see for example 
iKuzio de Nar av et al. 2006; de Blok et al. 2001) and giant low 
surface brightness galaxies (iPickering et al.lll997l) . On the con- 
trary, a much wider range of values is measured for regular 
galaxies, from few 10"^ Mq pc~^ to several 10^ M p pc"^ (i.e. 
ISalucci & Borrielldl200ll: iNoordermeer et al.ll2007h . Therefore, 
in this picture, the bulge of ESO 323-G064 resembles more 
the central mass density of regular bulges than those measured 
in low surface brightness galaxies. This is consistent also with 
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Fig. 10. confidence levels for the pseudo isothermal dark matter halo model. The crosses represent the location of the best fit in 
the parameter space (x^ = 1.5). 




Fig. 11. Fraction of dark matter mass compared to total mass, as 
a function of radial distance, for the inner 5" of ESO 323-G064. 
NFW (solid line); pseudo isothermal (dashed line). 



Fig. 12. Filled circles: mass density values derived from the ob- 
served velocity and velocity dispersion (Equation 12111 compared 
to the best fit NFW (continuous line), pseudo isothermal (dashed 
line) and self consistent (dotted line) model predictions. 



the fact that bulges of giant LSB galaxies are photometrically 
similar to those of regular high surface brigt htness galaxies 
jMcGaugh et al.lll995l;lBeiiersbergen etai][T999h) 

Acknowledgments. The authors wish to thank the referee A. 
Bosma for useful suggestions which improved the paper content 
and the discussion. 
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Appendix A: Constraining the halo parameters with 
disk i-iyS velocity curve and surface brightness 

The HyS velocity curve measured in the disk between 10" < R < 
30" (Section inT i and the disk surface photometry from the ESO- 
LV images can be used to constrain the dark halo parameters, 



under the assumption that the measured Hy6 velocity curve is a 
good representation of the galaxy circular velocity. 

A1. Circular velocity for the disk and the halo 

The circular velocity predicted for a pseudo isothermal halo is: 



Visoiir) = VhJI-— arctan — 



(22) 



where Vh - ^flAnGpor^) and po, '"o are the halo parameters de- 
fined in Section [573] 

The circular v elocity predicted for a NFW halo is 
(lNavarroetal.lll997h : 



yNFw{r) - Vaoo 



1 ln(l + cx) - cx/(l + cx) 
X ln(l + c) + c/(l + c) 



1/2 



(23) 



In Equation|23]c is the concentration parameter, related to the 
density parameter p, (used in our fit procedure) and the critical 
density pcrnic - 3H^/(^7tG) (H - 75 is the adopted value for the 
Hubble constant) with the following Equation; 



PsIPc. 



200 



3 ln(l H-c)-c/(l -Hc) 



(24) 



The parameter V200 in Equation |23] is the circular velocity at 
'"200 - {vs is the halo scale parameter used in our fit proce- 
dure) and it is defined as: 



V2()0 - 



/ GM2oo \ 
\ '"200 



1/2 



where 



An 



M20O = 200— ■ Pcrilwr200- 



(25) 



(26) 



Together with the halo circular velocities defined by 
Equations |22] and |23] we have to add the contribution given by 
the stellar disk (the bulge contribution is negligible at this radial 
range). Assuming an exponential disk with central mass surface 
density Mp and scale length h, the predicted circular velocity is 
(lFreemanlll970l) : 



Vdisk(r) = -v/47rGMo 



/o(|^)^o(|^)-A(|-)^,(|-) 



4/1 



(27) 



where /„ and K„ are modified Bessel functions of the first and 
second kind. 

To obtain the mass surface density Mq and the disk scale 
parameter we retrieved ESO-LV images of ESO 323-G064 in 
the B and R b ands and perform ed a photometric decomposition 
with galfit dPeng et al.ll2002l) . 

The result of the galfit decomposition are listed in Table 
|2] From the total disk magnitude in R and B band we compute 
a color o f (B - R)disk - 1-70 + 0.5. Using the prescription by 
[Bell & de JongI ( l200lh we can compute the mass-to-light ratio 
of the disk in the R band {MIL)disk = 4.81 + 0.77. The compu- 
tation made use of the relation between color a nd mass-to-light 
ratio for different models as listed in Table 3 of lBeif & de Jon^ 
( llOOlh . We use the relation in the R because the scatter between 
different models is smaller. 

From the central disk surface brightness //q = 21.2 ± 0.3 
mag arcsec"^ measured from the galfit decomposition in the 
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Table 2. Galfit best fit parameters. 





K-OQ.UQ 


D-Dana 


Bulge: 






Mtot 


lo.zJ ± U.J mag 


i /. / J ± (j.j mag 


n 

K, 


0.65 ± 1 arcsec 


0.55 ± arcsec 


b/a 


0.71 ±0.3 


0.71 ±0.3 


PA 


13 ± 20 deg 


-72 ± 20 deg 


Disk: 






M,o, 


14.11 ± 0.5 mag 


15.81 ± 0.5 mag 


Mo 


21.2 ± 0.4 mag 


23.3 ± 0.4 mag 


h 


12.4 ± 1 arcsec 


13.6 ± 2 arcsec 




0.53 ± 0.3 


0.52 ± 0.3 


PA 


33 ± 5 deg 


32 ± 6 deg 



Note: **' A disk axial ratio b/a = 0.53 corresponds to 
an inclination = 60° , assumin g an intrinsic disk axial 
ratio 01^0=0.18 jGuthriell993) . 

7?-band, assuming a distance of D - 194 Mpc, an extinction of 
Ar - 0.231 mag (from the NED database) the disk mass-to- 
light ratio (M/L)ijisk = 4.81 + 0.77 we calculate the central mass 
surface density Mq = (49 1+201 )M0 pc^^ = (4.3 + 1.8) ■ 10** Mq 
arcsec"^. 

Therefore, we can compute the total circular velocity Vc by 
adding the contribution of the light (Equation |27]| to the contri- 
bution of the halo (Equations|22]or|23]depending on the adopted 
scenario). 

A2. Mass model using stellar and gaseous kinematics 

We have also performed the mass model fit as described in 
Section|5]using the constraints from the H/3 circular velocity. As 
in the previous case we fit first the empirical bulge stellar rotation 
curve (to obtain rca and r/,, Equation [T2li and then we fit simul- 
taneously the stellar velocity dispersion (i.e the filled circles in 
Figure[8j and the HyS circular velocity (i.e the open diamonds in 
Figure |2]) data. 

Fit results are shown in Figure [13] and the best fit parameters 
are listed in Table |3] The results are similar to those in Table 
[1] With the additional constraints from the disk, both pseudo 
isothermal and NEW scenario give a reasonable fit to the data 
and it is not possible, given the errors, to say which is the best 
one. 

The formal errors given by the fit algorithm are slightly lower 
than the ones listed in Table [1] because of the addition of data 
points in the fit. However, the fit errors do not include the uncer- 
tainties from the disk photometry and the errors in the relation 
between mass-to-light ratio and disk (B - R) color. 

Errors in disk central surface brightness and in the mass-to- 
light ratio translate into an error of AMq/Mq ~ 41% on the mass 
central surface density, which is an eiTor of ~ 20% on the disk 
velocity curve {Vdisk °^ VM^)- However, the fact that the two 
different fit approaches lead to similar results is reassuring. 

Deep photometric observations and an accurate determina- 
tion of the mass-to-light ratio are highly desirable to better con- 
strain the mass distribution in ESO 323-G064. 



Table 3. Best fit parameters from the simultaneous fit of the stel- 
lar velocity dispersion and the H/3 rotation curve. Values are con- 
verted assuming D = 194 Mpc on the right column. 



Parameter Value 
Rotation curve: 

v„ 174 ±26 [km s-'] 

rn 1.6 ±0.8 [arcsec] (1.5 ± 0.8 [kpc]) 

Self consistent = 3-6 

Ml 6.6!|:« [10'OMo] 

NFW: = 0.9 

Ml O.54IO4+ [10'«Mo] 

log.op, 10.0!°:^ Pogio (Mo/"')] (p., = nt_l [Mo pc-3]) 

r, 0.85!° [arcsec] (0.80 ± 0.18 [kpc]) 

Isothermal: x^ = ^-^ 

Ml 3.5!^ [lO'^Mo] 

log.oPo 9.6!5!:2 [log,o (Mo/"3)] (po = 5!? [Mo pc-3]) 

ro 0.39:^5;;;^ [ai-csec] (0.37 ± 0.08 [kpc] ) 
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Fig. 13. Model results obtained including the disk data. Upper 
panel: stellar velocity dispersion (filled circles, see Section 
15.2b compared to the NFW {solid line), the pseudo isothermal 
{dashed line) and the self consistent {dotted line) best fit mod- 
els. Lower Panel: measured circular velocity {open circles, from 
the HyS rotation curve derived in Section 13.1b compared to the 
NFW {solid line), the pseudo isothermal {dashed line) and the 
self consistent {dotted line) best fit models. 



